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We demonstrate that semiclassical fluctuations, their relaxation, and the 

Oh! 

^ ■ chiral phase transition are automatically incorporated in the numerical sim- 
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ulations of the classical equations of motion in the linear cr-model when lon- 
gitudinal and transverse expansions are included. We find that domains of 
disoriented chiral condensate with 4-5 fm in size can form through a quench 
while an annealing leads to domains of smaller sizes. We also demonstrate 
that quenching cannot be achieved by relaxing a chirally symmetric system 

through expansion. 
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One of the proposed explanations for the Centauro events [|l| in high energy cosmic 
ray experiments is the coherent emission of pions from a large domain of disoriented chiral 
condensate (DCC) 0. However, if the system has to go through an equilibrium phase 
transition, quark masses, though much smaller than the intrinsic QCD scale, prevent DCC 
domains from reaching a size larger than 1/Tc [Q. As an alternative, Rajagopal and Wilczek 
01 proposed that a nonequilibrium phase transition through quenching can generate large 
DCC domains. Their numerical simulation indeed observed the amplification of the long 
wavelength modes of the pion fields, but did not shed light on the size of DCC domains. 
Similar simulations by Gavin, Gocksch, and Pisarski ||^ are, however, not conclusive about 
the exact size of DCC domains, since they only looked at the pion fields averaged over the 
transverse dimensions which cannot reveal domains smaller than the lattice size. 

To model a quench, Rajagopal and Wilczek argued that one can evolve the classical fields 
according to the zero temperature equations of motion from a chirally symmetrical initial 
condition with short correlation lengths. In this Letter, we shall argue that semiclassical 
fluctuations introduced in the initial configuration actually render the effective potential to 
a non-zero temperature one. The interaction between the mean fields and the semiclassical 
fluctuations as well as their evolution can be automatically included in the numerical simu- 
lations of the equations of motion. Using an ensemble averaging technique, we demonstrate 
the relaxation of the fluctuations and the occurrence of the chiral phase transition due to 
the longitudinal and transverse expansions which are consistently included in our study, 
whereas in earlier simulations the system only evolves in an approximately constant 
but non-zero temperature effective potential. By choosing different initial configurations 
(both the mean fields and the fluctuations), we study the evolution of the system in both 
quenching and annealing scenarios. 

In the standard linear cr-model, the equations of motion are given by, 

g^-V'U = Hv'-<P')<P + Hn^, (1) 

where (j) = (cr, tt) is a vector in internal space, n„ = (1,0), and Hrifj is an explicit chiral 
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symmetry breaking term due to finite quark masses. In the following, we shall use A = 19.97, 
V = 87.4 MeV, and H = (119 MeV)^ with which = 135 MeV and = 600 MeV. We 
emphasize that the (p fields in the above equations can include both the mean fields and the 
semiclassical fluctuations. Thus, we can separate into two parts, 

0=(0) + 50, (2) 

where (0) are the mean fields and S(f) are the semiclassical fluctuations around (0). Using 
Eq. d^) and taking the average of Eq. (P, we obtain the equations of motion for the mean 
fields in the Hartree approximation [H, 



- V^^ (0) = \{v' - (0)2 - 3(505) - + Hn.^ 



(3) 



where (0)^ = {(pi) {(pi) , 6(p\\ is the component of the fluctuation parallel to (0), and 6(p± is 
the orthogonal component. These equations imply that the motion of the mean fields is 
determined by an effective potential, 

vm) = ^m' + 3(505) + (^^i) - ^')' - ^(^)' (4) 

which clearly differs from the zero temperature one in the presence of the fluctuations. By 
varying the level of fluctuations, chiral symmetry can be restored or spontaneously broken. 
The above effective potential is very generic since no assumption has been made for the 
fluctuations except that they are only of a classical nature. If the fluctuation terms in Eq. 
(H) are replaced by their counterparts in a finite temperature field theory, the well-known 



one loop effective potential at finite temperature [|T0| is recovered. Therefore, for the mean 
fields, and as we shall show for the fluctuation fields as well, the classical equations of motion 
have already included the effect of fluctuations present in the effective potential. This might 
look surprising, but can be easily understood in a thermal equilibrium case. In a flnite 
temperature fleld theory, the temperature dependence of the effective potential arises from 
the on-shell part of the propagator. Since no contribution from virtual particles is involved, 
all thermal corrections at one loop level are purely classical. 



Since we shall neglect the corrections due to quantum fluctuations P,p!T|, we can also 



derive the equations of motion for the fluctuation fields 5(f) by subtracting Eq. from 
Eq. (|^). These two sets of coupled equations for (0) and 50 can, therefore, be consistently 
solved through numerical simulations of the classical equations of motion. Since the 
fluctuations and their squared ensemble averages can evolve with time according to a given 
relaxation mechanism |PJ^J^, the time evolution of the field configuration obtained from Eq. 
(H), already includes the effect of the time dependence of the effective potential. The use of 
equations of motion does not ensure that the effective potential takes its zero temperature 
form or that chiral symmetry is spontaneously broken. What matters most is the initial 
fluctuation of the system. 

When 5"^ = (3(50|) + is large enough, chiral symmetry is restored (approx- 

imately, due to if 7^ 0). If the explicit chiral symmetry breaking term is neglected, the 
phase transition takes place at the critical fluctuation, 5l = v"^ /Q. For 5"^ < 6"^, the effective 
potential takes its minimum value at (0) = {a^, 0), where depends on 5^. When the mean 
fields are displaced from this equilibrium point to the central lump of the "Mexican hat" 
{{(pi) ~ 0), modes below the critical momentum, 

X{v' - {^)^ - 3{S(Pl) - {S(Pl))f\ (5) 

become unstable and thus DCC domains can form. Since the domain size is directly related 
to the time scale during which these modes are unstable, it strongly depends on the initial 
condition, 6"^ and (0) of the system. 

Let us now consider three different scenarios, (i) In a quenching scenario, the initial 
fluctuation is below the critical value, 6"^ < 6"^, and {(pi) ~ 0. As the mean fields roll 
down from the central lump of the potential, pion modes below will be amplified. In the 
meantime, as the system cools down and the fluctuation decreases via, e.g., an expansion, 
the effective potential will also change and the equilibrium point of the potential (cre,0) 
moves towards the zero temperature value (/7r,0). This will increase the roll-down time 
and lead to a larger domain size, (ii) On the other hand, if we initially choose {(p) to be 



very close to the equilibrium point of the effective potential, (cre,0), the mean fields and 
the effective potential may both evolve so that the system always oscillates around the 
equilibrium position. We refer to this scenario as a cold annealing since the system starts 
with an effective potential in which chiral symmetry is spontaneously broken, (iii) What 
we shall call a hot annealing scenario is similar to the cold annealing except that the initial 
fluctuation is much larger than the critical value, 5^ ^ 5^, so that chiral symmetry is almost 
restored. In both annealing cases, the mean fields can evolve almost synchronously with 
the effective potential so that the system only oscillates around the equilibrium point, i.e., 
(0)^ ^ t>2 — 3(50|) — One can expect then that the low momentum modes are less 

amplified and the domain size is smaller than in the quenching case. 

In order to confirm our arguments, we have carried out numerical simulations of Eq. (|l]) 
including both the longitudinal and transverse expansions. We assume boost invariance in 
the longitudinal direction so that the longitudinal expansion is automatically included. To 
consider the transverse expansion, we use a cylindrical boundary condition. The initial (p 
fields are randomly distributed according to a Gaussian form with the following parameters: 









0, 


(-^) - W)' = 


(vrf) = 5oV(r), 


(^) = 


(vf.) = 0, 







(6) 

where r = (x^ + y^)^/^ is the radial coordinate, the dot stands for the derivative with respect 
to the proper time r = (t^ — z'^Y^'^, o"o and 5q are constants which we can vary for different 
scenarios. We have introduced an interpolation function 

'r-Ro 



fir) 



exp 



1 



(7) 



r 

to describe the boundary condition. Rq is the radius of the initially excited region where 
fluctuations exist and the mean fields are different from their vacuum expectation values. 
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Outside this region, the vacuum configuration, = (/^r, 0), is imposed. T is the thickness of 
the transient region. The results presented in this Letter are obtained with Rq = 5 im and 
r = 0.5 fm. 

In our numerical simulations, we have used two-step Lax Wendroff method, a version 
of the leap frog method |T^. We have used an initial correlation length, £corr = 0.5 fm. 



Usually the spacing a of the lattice on which numerical simulations are carried out has been 
identified to £corr MM- To include the initial correlation and to reduce the finite size effect, 
we have adopted a lattice spacing smaller than icon HHI- The initial fields are therefore 



uniform within icon x icon squares. Since domain formation is caused by the amplification 
of low momentum modes, the domain size should not be very sensitive to the value of icon- 
We have actually confirmed this in our numerical simulations [1^ . 



For the three scenarios we consider here, we take (i) 6^ = 6"^ = with which the 

system is about to go through a phase transition. The initial field configuration is set to 
ctq = for the quenching case, (ii) For a cold annealing case, we take cxq = = 44 MeV, in 
order to have an equilibrium initial configuration with the given fiuctuation. (iii) For a hot 
annealing case, we take 6^ = and (Jq = cTg = 20 MeV. 

We define a correlation function C (r, r) as 

where the sum is taken over those grid points i and j such that the distance between i 
and j is r. In Fig. 1(a), we show the time evolution of the correlation function for the 
quenching case. Throughout the calculations shown in this Letter, we take the initial time 
To = 1 fm and a lattice spacing a = 0.25 fm. At r = tq, there is no correlation beyond the 
initial correlation length £corr- We can clearly see that a long range correlation emerges at 
later times. Note that the apparent shrinking of correlation length at r = 7 fm is a little 
misleading. Since the resultant pion distributions only depend on nf, C{r, r) < should be 
regarded as the manifestation of the correlation as long as C(r, r) ^ 0. At r = 7 fm, the 
typical correlation length is as long as r ~ 2.5 fm. In Fig. 1(b), we compare the results of 
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the quenching, hot and cold anneahng scenarios at r = 7 fm. We observe that quenching 
gives the largest correlation among the three cases. We have also checked that changing 
o"o to in the hot annealing case does not help much to create larger domains, since the 
system has moved to the equilibrium position before the chiral phase transition takes place. 
In other words, a quenching condition can never be realized through a hot annealing. The 
situation does not improve much even if a second order phase transition is assumed {H = 0), 
because the expansion time scale is too short for any long range correlation to develop. In 
the quenching case, the expansion of the system reduces the fluctuation, and as a result, 
the evolving effective potential provides a longer roll-down time for the system to form 
larger domains. We have in fact checked that for larger tq cases, where relative longitudinal 



expansion is slower, a smaller correlation is generated [|1^ . 

In Fig. 2, we show the time evolution of the average a field, (a), and the average 
fluctuations, (5or2)V2 ^nd {Siriy^'^. {Sir^y/^ and (57r|)^/2 g^^g similar to (^vrf)^/^. In this 
calculation, the initial fluctuation, 6^ = 3f^/8, has been taken to be much larger than the 
critical value so that chiral symmetry is almost restored. For (Tq, we have taken its value 
to be (Te = 9 MeV. We have used 100 events and averaged over the central region r < 3 
fm. We see that (5(t^) and {Snf) decrease on the average with time due to the longitudinal 
and transverse expansion. On the other hand, (a) increases, following the equilibrium 
position of the evolving effective potential. In principle, (a) approaches to its vacuum value, 
as (5^ goes to zero. A very interesting and important point is that {Sa'^) and {STrf) 
decouple from each other at about 6 = 35 MeV, which is about the value of the critical 
fluctuation, 5c = {v'^/GY^'^. This decoupling is nothing but the manifestation of the mass 
splitting during the chiral phase transition: the pion mass becomes smaller and the sigma 
mass becomes larger. Obviously, the semiclassical fluctuations also experience an evolving 
effective potential as well as the mean fields ( as seen in Fig. 2, (a) rapidly increases after 
the phase transition as the system follows the ever-changing equilibrium position). 

Finally, to demonstrate domain formation, we show the contour plot of tt2 in Fig. 3 for 
one event in our quenching function of the transverse coordinates at the initial 



time r = To = 1 fm and r = 5 fm. We can clearly see two large domains with opposite signs 
in at r = 5 fm in this event, whereas initially there exists no structure inside the chirally 
restored region. We also note an apparent transverse expansion and decreasing fluctuations 
in the inner region. The domain formation is more dramatic for a smaller system . 

In summary, we have shown that the usual prescription for a quench actually already 
includes the effect of fluctuations. The relaxation of the semiclassical fluctuations and their 
effect on the effective potential are automatically included in the evolution of the system 
which undergoes both longitudinal and transverse expansions. Mass splitting of the pion 
and sigma flelds at the classical level during the chiral phase transition has been clearly 
demonstrated. We have also shown in our numerical calculations that DCC domains of 
a typical size up to 4-5 fm can form for realistic parameters in the linear cr-model if the 
quenching initial condition is realized. Furthermore, we have demonstrated that such a 
quenching condition cannot be achieved by relaxing a system from a chirally symmetrical 



phase through expansions. We will discuss elsewhere |T2[ whether and how a quenching 
initial state can be realized in hadronic or nuclear collisions. 
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Figure Captions 

Fig. 1 Correlation functions (a) at r = 1 fm (initial state), 3 fm, and 7 fm for the quenching 
initial condition; (b) for the quenching, hot and cold annealing scenario at r = 7 fm. 

Fig. 2 Evolution of (a), ((5(T^)^/^, and (^tt^)^/^ as functions of r. The initial fluctuation is 
taken to be, 5^ = Sv^/S, and ctq = 9 MeV. The averages are made over 100 events and 
within r < 3 fm. 

Fig. 3 Contour plot of 7r2 field in an event at r = tq = 1 fm and t — 5 fm. A quenching 
initial condition is used as in Fig. 1. The 7r2 field has opposite signs in the two large 
domains. 
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